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Abstract 

A general time-dependent projection technique is applied to the 
study of the dynamics of quantum correlations in a system consist- 
ing of interacting fermionic and bosonic subsystems, described by the 
Jaynes-Cummings Hamiltonian. The amplitude modulation of the 
Rabi oscillations which occur for a strong, coherent initial bosonic 
field is obtained from the spin intrinsic depolarization resulting from 
collisional corrections to the mean-field approximation. 

PACS number(s): 03.70+k, 42.50.Md, 32.80-t 
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1 Introduction 

Even very simple open subsystems of closed quantum mechanical systems 
can display very intricate dynamical behavior, described by an effective, non- 
unitary time evolution law. The quantum state of each such subsystem can in 
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fact be described in terms of a reduced density operator which will in general 
evolve non-unitarily on the account of correlations and decoherence effects 
involving different subsystems |TJ, |2], [tj . The non-unitary effects will manifest 
themselves specifically through the dynamical evolution of the eigenvalues 
of the subsystem reduced density matrices, so that individual subsystems 
evolve in general in a non-isoentropic fashion. 

A first mean-field-like approximation to this general picture consists in 
assuming isoentropic subsystem evolution under effective, time-dependent 
Hamiltonian operators for each subsystem. This approximation can be made 
self-consistent by allowing for explicit dependence of each subsystem effective 
Hamiltonian on the present state of the other, external subsystem [|TJ. Even 
when such selfconsistency requirement is abandoned by simply parametrizing 
the externally induced time dependence of the relevant subsystem effective 
Hamiltonian, one is still left in general with a nonlinear dynamical problem 
which defies exact treatment in a many-body or quantum field theoretical 
context. In this case, the subsystem time evolution can be analysed under a 
further mean field approximation through a variational principle by making 
a Gaussian Ansatz for the subsystem density matrix ||. 

In this work we implement the fully self-consistent double time-dependent 
mean field approximation to a simple model system consisting of interact- 
ing fermionic and bosonic subsystems described by the Jaynes-Cummings 
Hamiltonian []|]. Furthermore, by using the projection approach developed 
by Nemes and Toledo Piza [|IJ we obtain and evaluate explicitly corrections 
to this self consistent mean field description. These corrections appear in the 
form of suitable memory integrals added to the mean-field dynamical equa- 
tions. The resulting dynamical equations acquire then the structure of kinetic 
equations, with the memory integrals performing as collision terms which 
eliminate the isoentropic mean-field constraint. Relevance of the Jaynes- 
Cummings model in this context stems also from the fact that, besides being 
soluble, it can be seen as corresponding to the relativistic scalar plasma [H 
in zero spatial dimensions. The study of soluble models provides a clear un- 
derstanding of the physical phenomena involved and is useful in controlling 
various approximations necessary for the treatment of more realistic cases. 

The definition of the relevant variables and our approach for studying the 
dynamics of these variables are described in section 2. In section 3 we obtain 
the dynamical equations for these variables in the collisional approximation . 
Finally we show and discuss some numerical results, comparing to the exact 
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solution of the model. 



2 Gaussian Variables and Their Effective Dy- 
namics 

The Jaynes-Cummings Model (JCM) Q describes the simplest fully quan- 
tized version of a system consisting of two interacting quantum subsystems 
of different nature: a two level atom (fermionic system) and a quantized field 
mode (bosonic system). The corresponding exactly soluble Hamiltonian is 
given by 



where the annihilation and creation operators b and satisfy the boson 
commutation relation [b, tf] = 1 and the fermion operators a\ and a\ satisfy 
the anticommutation relations {a\,a\,} = 5\> t \. The solubility of the model 
stems from the fact that the interaction term of H is written in the so called 
Rotating Wave Approximation, which consists in ignoring "anti-ressonant" 
terms involving afa_i^ and its Hermitian conjugate. 

The general idea of our approach is to focus on the time evolution of 
variables which are fully determined by the mean values of the annihila- 
tion/creation operators and their bilinear or quadratic combinations, hence- 
forth refered to as gaussian observables. The state of the system is generally 
given in terms of a density matrix T in the Heisenberg picture. T is hermi- 
tian, time independent and has unit trace. 

The mean values of the boson operators are given by 



and its complex conjugate. Here and in what follows the symbol Tr BF denotes 
a trace over both bosonic and fermionic variables. Partial traces over bosonic 
or fermionic variables will be written as Tr B and Tr F respectively. Using them 
we can define the shifted boson operators 





(&(*)> =B t = Tr BF b{t)F 



d(t) EE b{t) - Bt 

which have vanishing ^-expectation values. 



(2) 
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Next, the mean values of the bilinear forms of boson operators can be 
combined in an extended one boson plus pairing density matrix 1Z |], 0] 



r b n B 
n* i + r e 



(dt(t)d(t)) (d(t)d{t)) 
(<fl(t)rf(t)) (d(t)rf(t)) 



(3) 



where the one-boson density matrix R B is hermitian and the pairing density 
matrix I1 B is symmetric. 

For the fermionic system the extended density matrix is 



rh 



-n; i-r f 



(a[,(t)a x (t)) (a y (t)a x (t)) 
(a[,(t)a\(t)) (a y (t)a{(t)) 



(4) 



where the hermitian matrix R F and the antisymmetric matrix I1 B are the 
one-fermion and pairing density matrices, respectively. 

To deal with the pairing densities one defines new Bogoliubov quasi- 
particle operators. In the case of the bosons we define quasi-boson operators 
as 

(3(t) = x* t (b-B t )+y* t (b + -B* t ) 

(5) 

/?+(£) = Xt (b + -Bl)+y t (b-B t ) 

and require that ((3(3) = ((3^(3^) = 0. The preservation of the commutation 
relations requires furthermore that the transformation coefficients Xt and yt 
be chosen so that |x t | 2 — \y t \ 2 = 1. Thus the densities (b'tf) and (bb) can be 
parametrized in terms of x t , y t and v = (ft (3). 

The coefficients Xt and yt of the Bogoliubov transformation are deter- 
mined by solving the secular problem 

GTZ B X B = X B GQ B (6) 

where 
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' xt y*t ' 




' 1 




' v 


x B = 


_ yt x* t _ 


G = 


.0-1 




_ 1 +v _ 



(7) 



The eigenvalue v can be interpreted as shifted boson occupation number for 
the paired natural orbital described by X B . Since the Bogoliubov transfor- 
mation is canonical one can verify that X B satisfies the orthogonality and 
completeness relation 



X^GX B — X B GX^ — G . 
For the fermion operators we define the transformation 



(8) 



ai(t) 
a\(t) 



X* 
X 



ai(t) 
a_i(t) 
a\(t) 



(9) 



where 



X 



u t -v t 



(10) 



Here we require that (ay(t)a\(t)) = (a\,(t)a\(t)) = and the fermion den- 
sities (a\ 

V\ = (a{ 
problem 



sities {a\,(t)a\(t)) and (ay(t)a\(t)) are parametrized in terms of u t , v t and 
Px = (a\a\), A = — 1, 1. As in the case of the bosons we have an eigenvalue 



where X v is the transformation defined in 



(11) 



X* 
X 



and Xl 



X T 

xt 



(12) 
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satisfying the unitary condition 



A'jA'p — fYpA'p — I4. 



(13) 



The matrix Q F is the extended density matrix in quasi-particle basis 



Pa 
1-Pa 



(al,(t)a x (t)) 







(a y (t)a{(t)) _ 



(14) 



where pa is the one-fermion density matrix in quasi-particle basis. In this 
basis Pa is diagonal with eigenvalues pi and p_i which correspond to the 
occupation numbers of the natural orbitals. 

The next step is to obtain the time evolution of the mean values of gaus- 
sian observables. We begin with the boson subsystem. For B t one finds 
immediately 



iBt = Tr BF [b, E\T = x t Tr BP [/3, E\T - y* t Tr BF [(3\ E\T 



(15) 



where H is the Hamiltonian given by d). For the remaining quantities, we 
first rewrite the eigenvalue equation (|), using (||), as 



X B 1Z B X B — Q i 



Taking the time derivative we have 



(16) 



xiii B x B = Q B -xln B x B -xln B x B 

= Q B — X^GX B GQ B — Q B GX] i GX B . 



(17) 



The left hand side of this equation can be evaluated using the Heisenberg 
equation of motion 



iX B TZ B X B 



I Tr BF [^(3,H]F Tr BF [PP, H]J- 
\Tr BF p{3\H]F Tr BF [(3(3\H]r 



(18) 
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The right hand side of eq. (|T7]), on the other hand can be rewritten using 
(0). Equating the result to (p~8|) we have equations which describe the time 
evolution of the bosonic variables explicitly : 



iv 



Tr^[fi(3,H]T 



(19) 



i{l + 2v){xlvt-x* t y* t ) 



Tr 



BF [HH 5 



H\T. 



(20) 



We have next to obtain analogous equations for the fermionic variables. 
Taking the time derivative of the eigenvalue equation ( |TT| ) and using the 
unitary condition (|13|) we obtain 



(21) 



Substituting (|1^) and ( |T4| ) and using the unitary condition (|T3"| ) we can eval- 
uate the right hand side of (|2l|). To evaluate the left hand side we use the 
Heisenberg equation of motion. We obtain 



IPX 



a\a x ,H 



-1,1 (22) 



i{Pi -P-i)(ikVt - v t u t ) 



ct{a:_i, H 



T 



(23) 



which describe the dynamics of the fermionic variables. 

Eqs. (JlJ), flT|), ©, (0) and (|23|) are however not closed, T is not fully 
determined by gaussian variables. In order to deal with this situation we 
split the full density as 



= ^Tl + T\t) (24) 

where the factorized form of !Fo(t) embodies what we refer to as the double 
mean field approximation. The subsystem densities JFq and JFq are in fact 
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gaussian densities, written in the the form of an exponential of a bilinear, 
Hermitian expression in the creation/annihilation parts of the bosonic and 
of the fermionic fields respectively. When explicitly normalized to unit trace 
they are given in terms of the transformed operators introduced in Eqs. (|^) 
and as § 



= n\p\a[a x + (1 - p x )a x a[] 
x 

1 / v \fftf> 
1 + u \l + u) ' 



•r o 



(25) 
(26) 



The remainder density T'it) is consequently a traceless, pure correlation part 
of the full density. In view of the special form used for JF (t) it will in general 
contain correlations of two types : inter-subsystem (boson-fermion) correla- 
tions and intra- sub system (boson-boson and fermion-fermion) correlations. 

The next step consists in using the fact that the decomposition ( |24|) can 
be implemented in terms of a time-dependent projection operator V(t) = 
V(t)V(t) (cf. Ref. §) 



V- = F*Tr F V B ■ +^Tr B V F ■ -^F*Tr F Tr E 



(27) 



where V B is the boson projection operator (11 



1 + 1/ J v{\ + V) V 



+ 



1 + u v ' 2u 
V F is the fermion projector 



2(1 + vY 



(28) 



TV 



A i^ 1:1 l-Px ) MAtli.A Pa/(1-Pa) J 

(29) 

and the dot stands for the object on which the operator acts. With the help 
of V(t) the factorized gaussian density can be obtained as 



•77F 
•'o 
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F (t) = V(t)F, 
while its time-derivative is given by the expression 



(30) 



& (t) = [F (t),H]+V(t)[H,F\ 



(31) 



which is in fact the Heisenberg Picture counterpart of the equation VT = 
which has been used to define V(t) in Schrodinger Picture ^ [|. 

The remaining formal steps towards closing the equations of motion for 
the gaussian variables are now straightforward. From fl24]) and (|3T| ) we can 
obtain a differential equation for T'it) in terms of J-'oit) 



F{t) = QCr Q (t) 



(32) 



where we introduced the complementary projector Q(i)J r = (X — V{t))!F = 
J-'(t) and the notation £■ for the Liouvillian time-displacement generator 
[H, ■ }. This equation has the formal solution 



(33) 



where the first term accounts for initial conditions, and Q(t,t') is the time- 
ordered Green's function 



G(t,t') = Texp \i [ drVMC 
. Jt> 



(34) 



Thus ^F'(t), and therefore J 7 , can be formally expressed in terms of ^(t 1 ) for 
t' <t and of initial correlations JF'(O). This finally allows one to express the 
dynamical equations (|15D , (|19D , (p0|), (p2|) , and (|23D as traces over functionals 
of J r o{t') and of the initial correlation part JF'(O). 



3 Approximation to the Collision Dynamics 

As can be seen from the formally exact Eq. (^), the correlation contributions 
to the dynamical equations involve traces over functionals of J-"o(1f) with 
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memory effects besides the initial correlations. Since these objects involve the 
complicated fully correlated time-dependence of the Heisenberg operators, 
approximations are clearly needed for their actual evaluation. 

A systematic expansion scheme for the memory effects in which the mean 
energy is conserved in all orders has been discussed in ref. |§ in terms of the 



Schrodinger Picture. Following Lin [11] we implement in the Heisenberg pic- 



ture a modified version of the lowest order approximation given in that work 
which consists in approximating the actual time evolution of the operators 
(3(t) when evaluating memory effects, by the effective mean field evolution 



= \j3, H ] + i{x* t x t - y* t Vt)f3 - i{x*y* t - x* t y* t )fi - Tr BF [f3, H]T. (35) 

The last three terms account for the explicit time dependence of the /3(t) 
related to the shift amplitudes Bt and to the Bogoliubov transformation (|5|). 
For the fermion operators a.\{t) and at_i(t) the corresponding approximation 
reads 



i ~d~t~ = \- ai,H °\ +^u t + VtV t )a 1 + i{u*v* t - u^)a^ 

(36) 

= [a-i, H ] + i(u t u* t + v t v* t )u-i - i(u t v t - u t v t )ai- 



In these expressions the Hamiltonian written as H (t) is taken as the effective 
mean-field Hamiltonian (cf. Ref. ||) 

Ho = Vi(t)H + flTr BF [(3,H]F'(t)-/3Tr BF [f3\H]F'(t) 

+ , a ' a ' 1 r Tr m [al iai , H]F\t) + ^ . Tr BP [a\^ H]T'{{$,7) 
(P-i-Pi) ' (Pi-P-i) 
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The first term in this expression correspond to the lowest approximation 
according to ref. ||. Note that V(t) is not an orthogonal projection, i.e. 
V'(t) 7^ V(t). The remaining terms, included here, represent correlation 
contributions to the effective mean-field. 

Consistently with this approximation, the Green's function fl34}) is in low- 
est order just the unit operator 

Q{t,t!)=Z 

so that the correlation density is approximated as 

T\t) = g(°\t,0)F'(0)-i f dt'Q(t')CF (t'). (38) 

Jo 

For simplicity, we restrict ourselves now to initially uncorrelated states, i.e. 
JF'(O) = 0. The dynamical equations of the gaussian variables with correla- 
tion terms then become 



[Tr BF {0 F {t)T)\ = Tr BF [O F (t),H]F (t)-iTr BF [O F (t),H) f dt'Q{t')CT (t' 

Jo 



d 

dt Jo 

mean-field approx. correlation term 

(39) 

4 [Tr BF (0 B {t)F)\ = Tr BF [O B {t),H]F {t)-iTr BF [O B {t),H] f dt' Q(t')£F Q {t') 
at Jo 

mean-field approx. correlation term 



where O f and O b correspond to fermion and boson gaussian variables, re- 
spectively. 

The terms involving To{t) can be calculated from H in quasi-particle basis 
and the projection operator V given by (PTQ, fl2"E| ) and fl2"9"|). However, the 
traces in the last terms of (|39|) still can not be taken directly, since the quasi- 
particle operators in the integral and in the first commutator are at different 
times. Therefore, we use the approximation in which the time evolution of 
these operators is given by ( |35l) and (|^). Using ([37]), (|2~0"1) and ( [231) we obtain 



.dp 
'~dt 



i(x* t x t - yly t )l3 + u{\x t \ 2 + \y t \ 2 )p = F (t)P 
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&0L\ 
I—— 

dt 

.dot— i 
dt 



-(\u t \ 2 -\v t \ 2 )-\(u* t v t B t + u t v* t B* t 



-(\u t \ 2 -\v t \ 2 )-\(u* t v t B t + u t v* t B* t ) 



at + i{u* t u t + vlv t )oti = F ai (t)ai (41) 

a_! + i(titu* + i) t Vf)a-i = F Q ,_ 1 (i)a_i 

(42) 



From (0), ([4l ) and ([42 ) one can verify that the operators /3(t), a(t) at 
different times are related as 



(43) 



where the phase (p(t,t f ) is given by 



<p(t,f) = - drF p {r) 



and 



z— = F Ql (t)ai 



Ql (t) = e ^(*'*') ai (t') 



a_ x {t) = e* a W> a_i(t) 



(44) 
(45) 



where the phase £ a (t, t') is 



= - [ drF a (r). 
Jt' 



Finally, the equations of motion that describe the present approximation 
to the collisional dynamics are 



v 



X 2 {(u* t v t y* t - utv*x* t )h + (ufy* t + vfx* t )h + (u 2 x* t + v 2 y*)I 3 ] + 



+ C.C. 



(46) 
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^(Pi-P-i) = 2\ 2 {-(ufx t + vfy t )i; + (ufy* t +v?x* t )I 2 } 

+ c.c. 

i(jpt - p-i)(u t v t - Utv t ) = (pi - p-i)[eu t v t + X(u 2 t B* t - v 2 B t )] + 

- 2i\ 2 (u* t v t x t - u t v*y t )I^ - 2i\ 2 (u t v* t x* t - v*v t y*)I 2 + 

+ i\ 2 (u 2 y t + v 2 x t )I* + i\ 2 (u 2 t x* t + v 2 y* t )h 

i(l + 2v)(x* t y* t -x* t y* t ) = -(1 + 2v)2ux* t y* t - 2i\ 2 (u* t v t y* t - u t v* t x* t )I* + 

+ 2 l \ 2 {u?y; + vfx* t )P 3 + 2i\ 2 {u 2 t x\ + v 2 y* t )I* 

where 

h = f dt'{u t ,v* t ,y t , - u;,v t ,x tl )\pi{l - Pl ) +p_i(l -p-i)]^ 1 '^ (50) 
Jo 

Jo 

h = f dt\ufx tl + vfy tl )\p x {l - p_ x ) + (pj - p_ 1 )^]^ i(2 « ( ''*' ) -^52) 



o 



Using the phase equation (|43"D one can write the integral I\ as 

Jj = e - I "j( 1 )(t) (53) 



where 



Ip(t) = I dt'Fpit') (54) 
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I«(t) = / ^'(^<,^-<,^x t 0[pi(l-pi)+P-i(l-p-i)] t 'e^'^5) 



i(a£x t - ylvt) + ^(|x t | 2 + \y t \ 2 ) 



(56) 



To solve the equation (|53| ) we write differential equation for the integrals 
1/3 and IP 



ip = i{x* t x t - y* t y t ) + u(\x t \ 2 + \y t \ 2 ) 



(57) 



/(!) = (w^-w^^fp^l-p!) +p_ 1 (l-p_ 1 )] t e- J ^ (58) 



Using the other phase equation we obtain, from the integrals I 2 and J 3 , 
the following differential equations 



j e (2) = KV + ^)[p-i(l-pi) + (p-i-pi)^] t e- i{2/Q+7 ' 3) (59) 



4 (3) = + vfy t )\pi(l - P-i) + (pi - P-i)u] t e 



(60) 



i(u*u t + v*v t ) + -{\u t \ 2 - \v t \ 2 ) - \(u*v t B t + u t vlB* t ) (61) 



Thus, we can integrate the equations (|57|-|5"S|) and (|59| - |6TD together with 
the equations ( ^6| - |49|) . These equations describe the time evolution of the 
gaussian variables in the collisional approximation. 



4 Results and Discussion 

We give below results obtained for two relevant observables associated with 
the spin subsystem of the Jaynes-Cummings Hamiltonian: the atomic in- 
version, (cr 3 ) t , and the "intrinsic" atomic inversion, defined in terms of the 
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projection of a along a unit vector pointing along the direction of the spin 
vector, (a)/\(a)\, which will be denoted by (cr p )t [0- The latter quantity cor- 
responds in fact to the degree of polarization of the spin 1/2 subsystem, and 
hence is directly related to the purity of the corresponding reduced density. 
Thus any time dependence of this quantity implies non-isoentropic behav- 
ior of the subsystem dynamics. In terms of the Bogoliubov parameters, the 
mean values of these observables are given by 

(0-3)* = (pi -p-i)cos29 

(<7p)t = (pi-P-l) 

where we use the parametrization 



u t = cos 9 



v t = e * s sm9. 

A typical example of the results obtained within the present approxima- 
tions for the time evolution of these objects is shown in the figures below. 
They have been obtained by solving numerically the dynamical equations (|46f - 
f49|) and are compared with the corresponding exact solutions of the model. 
Figure 1 shows the time evolution of the atomic inversion (03) t when the 
atom is initially prepared in the excited state (pi(0) = 1, p_i(0) = 0) and 
the field is initially in the coherent state constructed as the vacuum of the 
displaced annihilation operator b (x t =o = 1, Vt=o = and z/(0) = 0). The 
mean photon number B^ =0 is 25 in this example. We use units such that 
e — oj — 1 and A = 0.5 . The evolution of the intrinsic polarization for the 
same initial conditions is shown in Figure 2. 

The mean-field approximation amounts to ignoring the integrals I±, I2 and 
^3, which are set to zero at all times in this case. This causes the occupation 
probabilities z/, pi and p-\ to become time- independent, so that there are no 
correlations developing between the subsystems. As a consequence, there is 
no depolarization and the inversion undergoes Rabi oscillations with constant 
amplitude. In the collisional approximation a modulation of the amplitude 
of the Rabi oscillations is obtained from the depolarization resulting from the 
approximate collisional corrections to the mean field approximation, included 
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in Eqs. ( |46| ) to (|52|), which allow for intersubsystem correlations to develop. 
The effectiveness of the present treatment of these corrections can be judged 
by comparing its results with the exact solution. As can be seen in the figures, 
the approximate solution is quantitatively adequate up to times comparable 
with the depolarization time. 

We may infer from this that the inclusion of the present approximation of 
the collision integrals in the dynamical equations is not only fundamental to 
generate the qualitative behavior associated with the decorrelation process 
related to the initial damping of the Rabi oscillations, but also successfully 
describes such effects quantitatively over a time span covering at least several 
Rabbi periods. Improvements of the collisional approximation implemented 
here are however needed if one wishes to extend time range of quantitative 
reliability. 

Finally, we briefly comment on the extension of the present treatment 
to 3+1 dimensions in a field-theoretical context (the scalar plasma ||). In 
general, an entirely analogous procedure will apply also to this case, if only 
one expands the field operators in terms of the eigenf unctions of the corre- 
sponding extended one-fermion and one-boson densities. This is in general 
a time-dependent basis whose evolution will be determined by generalized 
forms of Eqs. ( f48[) and (f49"l). The case of a spatially uniform system is par- 
ticularly simple, since in this case the extended densities are diagonal in a 
plane-wave representation. It can be handled therefore in complete analogy 
with the 0+1 dimensional case if one uses a momentum expansion of the field 



operators [fn]| . As expected, the usual infinities of quantum field theory show 
up and can be absorbed by introducing appropriate counterterms. Work in 
this line is under way. 
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Figure Captions 



Fig.l. Time evolution of the atomic inversion with A = 0.5. Initial 
conditions: v — 0, x t — 1, fjt — 0, \B t \ = 5 (coherent state); pi = 1, p_i = 0. 
Full line: exact solution; dot-dashed line: collisional approximation; dashed 
line: mean-field approximation. 

Fig. 2. Time evolution of the intrinsic polarization. The initial conditions 
are the same as Fig. 1. Full line: exact solution; dot-dashed line: collisional 
approximation. 
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Figure 1 




Figure 2 




